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Abstract 

We express the heavy quark diffusion coefficient as the temporal variation of a Wilson line along 
the Schwinger-Keldysh contour. This generalizes the classical formula for diffusion as a force-force 
correlator to a non-abelian theory. We use this formula to compute the diffusion coefficient in 
strongly coupled M = 4 Yang- Mills by studying the fluctuations of a string in AdS§ x S$. The 
string solution spans the full Kruskal plane and gives access to contour correlations. The diffusion 
coefficient is D = 2/y/XTrT and is therefore parametrically smaller than momentum diffusion, 
7//(e + p) = 1/AttT. The quark mass must be much greater than T\f\ in order to treat the quark 
as a heavy quasi-particle. The result is discussed in the context of the RHIC experiments. 
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I. INTRODUCTION 



The experimental relativistic heavy ion program has produced a variety of evidences 
which suggest that a Quark Gluon Plasma (QGP) has been formed at the Relativistic 
Heavy Ion Collider (RHIC) 0, 0- One of the most exciting observables is the medium 
modifications of heavy quarks. In particular the electron spectrum from the semi-leptonic 
decays of heavy quarks is substantially suppressed relative to scaled proton-proton collisions 
0, 0| . Furthermore preliminary measurements indicate that the heavy quark elliptic flow is 
significant although less than the light hadron elliptic flow 0]. 

A variety of phenomenological models have estimated how the transport mean free path 
of heavy quarks in the medium is ultimately reflected in the suppression factor and elliptic 
flow 0, 0, 0|. The result of these model studies is best expressed in terms of the heavy 
quark diffusion coefficient. (In a relaxation time approximation the diffusion coefficient is 
related to the equilibration time, r^ eavy = ^D.) There is a sense from the models that if 
the diffusion coefficient of the heavy quark is greater than 

D> 1 - 

^ ~ rp > 

the heavy quark medium modifications will be small and probably in contradiction with 
current data. This interpretation of the RHIC results is perhaps too naive since the diffusion 
coefficient dictates the dynamics of non-relativistic heavy quarks. Diffusion may be irrelevant 
for the dynamics of the mildly relativistic heavy quarks measured at RHIC where radiative 
energy loss may be significant (UEll- Nevertheless, the diffusion coefficient is a fundamental 
parameter of the plasma and is essential to any discussion of the RHIC heavy flavor data. 

Unfortunately, for the experimentally relevant range of energy densities the QGP is not 
weakly coupled, and it is not easy to determine this transport coefficient. Ideally, the dif- 
fusion coefficient should be measured on the lattice, but this is difficult 11, 12f • Most 



theoretical works either compute the diffusion coefficient in perturbation theory and sub- 
sequently extrapolate to strong coupling (see e.g. 0), or develop models for the strongly 
interacting QGP |l3[- In this work we will compute the diffusion coefficient in strongly 
coupled M = 4 Super Yang Mills (SYM) where rigorous computations are possible if the 
AdS/CFT conjecture is accepted. Although this theory is not QCD, computations in this 
theory serve as a foil to the extrapolations based on weak coupling. 

This conjecture states that for a large number of colors, strongly coupled M = 4 SYM is 
dual to classical type IIB supergravity on an AdS^ x S 5 background 14 , H, 16||. (For reviews 
and lectures see Re f. (ljl. laP T he p hysics of heavy quarks has been studied using semi- 
classical strings On the gauge theory side a heavy charge is realized 
by breaking the U(N) gauge theory to U(N — 1) x U(l) through the Higgs Mechanism. 
The resulting W boson transforms in the fundamental and is heavy if the scalar expectation 
value is large. On the gravity side this corresponds to placing one of the D3 branes far from 
the remaining N — 1 branes. The dynamics of the heavy quark is dictated by the classical 
dynamics of the Nambu-Goto string stretching in the AdS$ x 5*5 background. The first thing 
computed was the expectation value of the Wilson loop to find the heavy quark potential 
(ljf . The result was immediately extended to finite temperature 21, 2^, 23[ 24 1. Since then 



there have been numerous studies of the other properties of Wilson Loops both at zero and 
non-zero temperature 17 . 
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One of the most interesting transport properties of the Af = 4 plasma is the heavy 
quark diffusion coefficient. Various transport properties have been computed using the cor- 
respondence leading to the remarkable conjecture that the shear viscosity to entropy ratio is 
bounded from below by 1/Att (111, EH Hi > I§t • This bound is influencing the interpretation of 
heavy ion results. Generally these calculations of transport in Af = 4 use the supergravity 
approximation. In contrast the computation of heavy quark diffusion will utilize the classical 
string theory directly. The diffusion computation in Af = 4 SYM is therefore experimen- 
tally and theoretically significant. For other applications of the AdS/CFT correspondence 



motivated by RHIC physics see Refs. |29|, |3C 

The computation proceeds as follows. In Section |H] we review the Langevin dynamics 
of the heavy quark. In Section IIHI we show that the strength of the noise in the Langevin 
process is determined by real time electric field correlators once the heavy quark has been 
integrated out. These electric field correlators can be expressed as fluctuations of a Wilson 
line running along the Schwinger-Keldysh contour. These results generalize the classical 
formula for the diffusion as a force- force correlator to non-abelian theories. In Section ITVl we 
compute the fluctuations of the string configuration which corresponds to this Wilson line 
running along the contour; the relevant string spans the Kruskal plane. Finally in Section IVl 
we summarize our results and discuss the implications for the RHIC experiments. 



II. LANGEVIN DYNAMICS 

In this section we will discuss the predictions of the Langevin process for determining the 
retarded correlator. We will assume that the Langevin equations provide a good macroscopic 
description of the thermalization of heavy particles 

d/Xi Pi 

~dt ~ ~M ' 
dpt 

dt 



The drag and fluctuation coefficients are related by the Einstein relation 

no = ^ . (2.2) 

The drag coefficient rjD can be related in turn to the diffusion coefficient 

T 2T 2 

D = — = — . (2.3) 
Mr] D k 



For a brief review of the Langevin equation and a derivation of these results, see Ref. |6[. 

The relevant time scale for medium correlations in a relativistic plasma is ~ D which is 
short compared to the relaxation time of the heavy particle (M/T) D. Furthermore, over the 
time scale of medium correlations the quark moves a negligible distance, v/T/M D. Thus, 
for the purposes of calibrating the noise (k), the mass may be taken to infinity and the 
heavy quark may be considered fixed. 
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More generally, the microscopic equations of motion for a heavy particle in the medium 

are 



f = ^(t). (2.5) 

Compare the response of the Langevin process to the microscopic theory. Over a time 
which is long compared to medium correlations D, but short compared to the time scale of 
equilibration (M/T)D, the drag term may be dropped, and we find 

J dt J dt' (&(t)fj(0) = ( time ) x kS H = j dt J dt> • ( 2 - 6 ) 
Taking the force in the ^/-direction for instance, we have 

dt (F y (t)F y (0)} . (2.7) 



We will drop the "y" in what follows. In QED this analysis relates the Langevin noise to an 



an electric field correlator |31|, |32| . The next section generalizes this result to gauge theories. 



III. MOMENTUM DIFFUSION IN GAUGE THEORIES 

We study the propagation of a heavy quark in a thermal bath of SU(N) gauge fields. After 
a Foldy-Wouthuysen transformation that reduces the heavy quark field to a two dimensional 
spinor field Q [3J|, the heavy quark effective Lagrangian is 

C = Q ] (id t - M -Aq)Q. (3.1) 

From this Lagrangian it is natural to identify the force in the Langevin equation with the 
operator 

T = J d 3 xQ\t,x)T a Q(t,x)E a (t,x), (3.2) 

where E is the chromo-electric field in the y direction. The momentum diffusion coefficient 
k is given by Eq. ()2.7|) 

k = J dt(F(t)f(0)) HQ , (3.3) 

where the average should be understood as a thermal average in the presence of a heavy 
quark - the meaning of this average is clarified below and in Appendix |X] 

In the real time formalism of thermal field theory (see e.g. EH)> the quantization 
of the fields is performed by extending the time coordinates into the complex plane along 
the Schwinger-Keldysh contour - see Fig. ^ There are two types of quantum fields: type 1 
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FIG. 1: The Schwinger-Keldysh contour (with T — > oo). Fields evaluated along the real axes 
are labeled as type 1, while fields evaluated on the —i(3/2 axis are labeled as type 2. The crosses 
indicate the insertion points of the color singlet force operator T (Eq. 13.2(1 . After integrating out 
the heavy quark, the crosses indicate insertions of the electric field, and the dotted lines indicate 
the path of the corresponding links. The electric field insertion may be rewritten as a variation (at 
times and t) of the Wilson line running along the whole contour. 

fields, evaluated on the real time axis and type 2 fields, which are evaluated on the second 
— i/3/2 time axis. We therefore define the correlation functions 



Temporal correlators of the operator T in the presence of an external heavy quark verify 
the Kubo-Martin-Schwinger (KMS) relations, provided single particle states may be consid- 
ered complete, i.e. gauge fields are not so strong as to create heavy quark-antiquark pairs 
and the density of heavy quarks is small - see Appendix [XJ In particular we define the 
Retarded Green function 




(3.4) 
(3.5) 
(3.6) 

(3.7) 



where T and T denote time and anti-time ordering respectively, and 




(3.8) 
(3.9) 



iG R (t)=0(t)([F{t),r{0)]) 



HQ i 



(3.10) 



and relate the remaining Green functions by inserting complete sets of states 

iG n (w) = iReG R (uj) - coth (|^) 1mG R (u) , (3.11) 

iG 12 {u) = iG 21 {u) = - — haG? a (o;), (3.12) 

iG 22 M = -iRe - coth (|^) ImG^u;) . (3.13) 

Using the fluctuation dissipation theorem, k is related to zero frequency limit of the retarded 
correlator 

2T 

re = lim ImG fl (cj) . (3.14) 

u>— >o o; 

This can be determined from iGi 2 (u;) for instance. 

For a Yang-Mills theory with external heavy quark the partition function is 



'HQ 



det [id^ — M — A — ie) 
x [[DQ^DQ^fc^i^c-M-^QQ^-iftQl^-T)) , (3.15) 

J I YM 

where the integration is performed along the Schwinger-Keldysh contour. 

The heavy fermion may be integrated out of this expression of the partition function. The 
propagator of the heavy quark in a fixed gauge background is computed from the Green's 
function 

i-^- - M - A, - ie^j G (x, t c] x', t' c ) = i6 (x - x') 5 (t c - t' c ) . (3.16) 

With this propagator, the partition function of the heavy quark yields the Polyakov loop 
along the contour 1 

Z HQ = V ps e-? M (W c [0]) ■ (3.18) 

Usually the contour of the heavy quark partition function is taken straight down the imag- 
inary axis to —i/3, yielding the usual Polyakov loop. However this is unnecessary and the 
answers are the same. 

Similarly, the heavy fermion fields may be integrated out of contour ordered force-force 
correlators 

(Ta[r(t a m0)]) BQ = (tr[U(-T - t{3, t c ) E(t c ) U(t c , 0) E(0) U(0, —T)]) , (3.19) 

where the tranporters U(tf,U) are calculated along the Schwinger-Keldysh contour with 
and tf the initial and final contour times. The denominator (Wc[0]) stems from the partition 



1 V ps is the phase space volume 

V PS = vJ-^ = V6*(0). (3.17) 
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function Eq. (|3.18|) . This correlator of electric fields can be generated by starting with a 
Wilson loop running along the contour and inserting small notches at contour times t — 
and t = tc as shown in Fig. [TJ 

Generalizing this procedure, we first introduce a source 5y(tc) for the force operator 

1 Z HQ \5y} = (T c e^o^)nt)\ . (3 . 20 ) 



Zhq[0) ^ l " J \ ° /HQ 

This is accounted for by modifying the heavy quark Lagrangian 

C HQ = Q ] {id tc -M-A - 8y{t c ) E{t c )) Q . (3.21) 
As in Eq. ()3.19|) . we integrate the partition function in the presence of a heavy quark 

Z HQ [6y] = V ps e-? M (W c [Sy}) , (3.22) 

where 

W c [Sy] = T c exp j-i jf dt (A + 5y(t)E(t))j . (3.23) 



This object can be easily expressed as 2 

W c [5y] = T c exp{-i J dtA ^T c exp^i J dt Sy(t)E(t)j , (3.24) 

where E is the dressed field strength 

E(t c ) = U(-T, t c ) E(t c ) U(t c , ~T) , (3.25) 

i.e. the transporter starts at an initial time — T, runs along the Schwinger-Keldysh contour 
from the initial point up to the contour time in, and returns to the initial time, — T. By 
means of the non Abelian Stokes theorem |3j| or some thought, it can be shown that 
Eq. ()3.24j) is the Wilson line along the time contour with a deformation in the y direction 
given by 5y(t c ) 

W c [5y] = T c exp j -i j dt {A + 5yA y ) J . (3.26) 

In summary, the path 5y(t) of a Wilson loop running around the Schwinger-Keldysh contour 
is the source for contour ordered force operators, T. In real time thermal field theory 
it is custumary to break up the source into 1 type sources and 2 type sources with the 
understanding that variations of the vertical part are not considered. The source Syi(t) and 
5y2(t) are variations of the the Wilson loop on the 1 and 2 axis respectively 

^ re i/*^(t)^(t) fe-^ dt ' s ^^) HQ = 1 <Wc[SViM) ■ (3.27) 

The momentum diffusion coefficient may then be written 

^ 1 / PW c [5yi, Syt] 



k = lim / dte +tuJt F w ( - ,\'\ ) ■ (3-28) 

<^oJ (Wc[0,0]) \ 6y 2 (t)6 yi (0) 1 



2 Sy is assumed to be small. 
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A. J\f = 4 Super Yang Mills 



The previous discussion was performed in a SU(N) gauge theory. However, we are in- 
terested on studying the interactions of a heavy W boson in a Af = 4 U(N) gauge theory 
broken to U(N — 1) x U(l). The heavy W boson propagator is given by the following Wilson 
line Q 

U = T c exp | -i J ds (A^ + \x\ &X 1 ) J , (3.29) 

where in addition to the gauge fields the adjoint scalars X 1 couple to the W boson through 
a Yukawa type interaction. Here G 7 is the VEV angle. 

The extension to M = 4 SYM amounts to changing the transporters and Wilson lines in 
accord with Eq. 13.291 In particular, using the non-Abelian Stokes theorem the M = 4 
analog of Eq. (j3~2^jl 

W c [Sy] = T c exp j-z J dt (A + O'X 1 ) J T c exp |z J dt 5y{t) £ S Y M (t) J . (3.30) 



Here we have defined 



Ssym = E + D t (Q 1 X 1 ) - D y (Q 1 X 1 ) , (3.31) 



and dressed it with transporters as in Eq. (j3.25|) . The necessity of the scalar derivatives 
is obvious after computing a plaquette with the links defined in Eq. ()3.29j) . In the low 
frequency limit the time derivatives of the scalars can be dropped and the QCD-like field 
strength is modified by gradients of the scalar fields. 



IV. THE ADS/CFT CORRESPONDENCE 

Fundamental Wilson Loops were first computed by Maldacena whose work we recall . 
On the gauge theory side heavy gauge bosons are introduced through the Higgs Mechanism. 
Specifically the U(N) gauge group is broken to U(N — 1) x £7(1) by giving a large expec- 
tation value to one of the scalars in the theory. The resulting W boson transforms in the 
fundamental representation of the remaining U(N — 1) gauge group. The equation of motion 
for the W boson is 

(id t - M - A - 6 1 X 1 ) W = , (4.1) 

and naturally leads to the Wilson loops considered in the previous section. 

On the gravity side this Higgs construction corresponds to placing one of the D3 branes 
far from the remaining N — 1 D3 branes. The propagation of the heavy W bosons is 
represented by semi-classical strings which are governed by the Nambu-Goto action 

S = — i- / drda^-detGuN d a X M d b X N . (4.2) 
zna J 

In the original Maldacena computation, the two parallel Wilson lines corresponded to a 
string world sheet that circumscribed the trajectory of a u>(x) and u>(y) bosons on the 
boundary. 
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FIG. 2: Kruskal diagram for the AdS black hole. The coordinates (t,r) span the right (R) 
quadrant. The dotted lines and the dashed hyperbolas represent the future and past horizons and 
the singularities, respectively. The thick hyperbolas on the sides of the two quadrants are the 
boundaries at r = oo. The Wilson line running along the Schwinger-Keldysh contour runs along 
along the "1" axis (the R boundary) and "2" axis (the L boundary). This corresponds to a string 
whose endpoints follow these boundaries. The minimal surface with these boundary conditions is 
the full Kruskal plane. 



Finite temperature is introduced by inserting a black hole into the AdS space. When 
the black hole is viewed in Kruskal coordinates there are two boundaries which are asymp- 
totically AdS. These boundaries correspond to the type "1" and "2" axis of the real time 
partition function (jjiL . Naturally the thermal description of a broken gauge theory at 
finite temperature also has w bosons running along the "1" and "2" axis. 

Since we are considering a single heavy gauge boson which propagates along the full 
Schwinger-Keldysh contour, there is a single Wilson line on the "1" axis and a single Wilson 
line on the "2" axis. It is natural to look for string solutions in the full Kruskal plane. Such 
solutions connect the type "1" w boson at the boundary of the right quadrant with the type 
"2" boson at the boundary of the left quadrant - see Fig. El 

To start, we look for string solutions in the first quadrant. The metric is 

ds 2 = ^ [-f(r)dt 2 + eJxf] + dr 2 + R 2 dQ 2 5 , (4.3) 

where f(r) = 1 — (y) 4 , R is the AdS radius, and r Q is related to the Hawking temperature 
irTR 2 = r . We define the scaled units, nTt = i, nTx = x and f = r/r and the metric 
reads 

^ = -f 2 /( f) dt 2 + f 2 d^ + 1 0- 2+d nl (4.4) 
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In what follows we will drop the "bar" until Eq. (|4.19|) 
Parameterizing the string as 

(r,cr) i-> (t = r, X||(r,(x), r = a, Q 5 (r,cr)) , 

It is straight forward to show and well known that 

X|| = Const and ^5 = Const , 

is a solution to the equations of motion. This solution extends to the full Kruskal plane and 
is the string solution which we are looking for. 

To see this we first write the metric in Kruskal coordinates 3 

ds 2 

— = g k (-dt 2 k + dx\) + r 2 dxjj + dQ 2 5 , (4.5) 

where the Kruskal scale factor is g k = (1/4) r 2 f(r) e~ 4r *. Parametrizing the string as 

(r,a) i-> (t k = r, X||(r,cr), x k = a, fi 5 (r,cr)) , (4.6) 
we determine the induced metric 

C = ^Z^th = R 2 [g 2 + g k (r 2 (xj,) 2 - r 2 (x,|) 2 + (fi' 5 ) 2 - (fi 5 ) 2 )] ^ , (4.7) 
where, for instance, xy = d T x.\\ and xj, = c^xy. Then we compute the canonical momenta 

_5C_ = g k r 2 ± ll 

(4.9) 



8C 9k r 2 x(. 



5(x[) v^ditA' 

and we see that the equations of motion 

„8C 5C . 
o T — + d a — =0, 

OX 0Xj| 

are satisfied trivially by X11 = Const. The canonical momenta vanish. A similar remark 
applies to the coordinates. Thus 

xii (tfc, rcjfc) = Const, VL^{t k) x k ) = Const, (4-10) 

is a solution to the string equations of motion throughout the Kruskal plane. This represents 
a Wilson line along the "1" time axis and back along the "2" axis of the Schwinger-Keldysh 
contour. 



Our notation closely follows Ref. |4|j], appendix A. Define r* = JJ" ^fhy\ + if so that r* is real in the 
right quadrant. The Kruskal coordinates are — e 4r * = tk~ x k and 2i = tanh -1 (t k /x k ) 
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Having constructed the string solution representing a W boson propagating all along the 
Schwinger-Keldysh contour we will proceed to fluctuate the shape of this string solution. 
Following the general philosophy of the AdS/CFT correspondence we equate the exponential 
of the classical action of the source to the generating functional 

_ e iS NG [Syi ,Sy 2 ] _ \ (W\8vi 6vo\) (A 11) 

e iS NG [ofi] e (W[0,0]) { [ VU m) ' { ' 

Now that the problem has been formulated along the Schwinger-Keldysh contour and KMS 
is satisfied, the results of Son and Herzog |38| may be adopted mutatis mutandis. (For 



related work, see Ref. |3flL l4lL l42j|.) For the retarded correlator 



POO 

G R (u) = -i / dte + ^6(t) {[F(t),r(p)]) HQ , ( 4 - 12 ) 
Jo 

the procedure is the following. In the right Kruskal quadrant, solve the classical equa- 
tions of motion for the source 5y(t, u) with infalling boundary conditions at the horizon. 
The retarded correlator may be found by evaluating the classical action of the string and 
differentiating twice with respect to the the boundary source, 5y(t,u = 0). 

Returning to the first quadrant with t,x\\,u= \ coordinates, the relevant string action 
with small fluctuations in the y-direction is 



R 2 r dt du 
NG ~ M / 2u^ 



i (yf\ 



l ~2 (j-^Wl)' 



(4.13) 



Notice that the infinite action corresponding to the unperturbed string is to be subtracted 
since it appears in the numerator and denominator of Eq. (|4.11|) . Fluctuating in the y 
direction, we define 

y(t,u) = J e- imt y(to)Y w ( U )^, (4.14) 
and the equation of motion of the fluctuating string is 



(2 + 6u 2 ) n tv 2 
dlY w - L 1 ^5 tt F ro + —Y n = . (4.15) 



This equation has regular singular points at the boundary u = and at the horizon 
u = 1. Near the horizon the solution behaves as 

Y m (u) = (1 - u) ±in '\ (4.16) 

Remembering that u = 1/V 2 , these solutions correspond to outgoing and infalling fluctua- 
tions. For the retarded propagator, 

/■oo 

G R (u) = -i \ dte + ^6(t) ([F(t),F{0)]) HQ , (4.17) 
Jo 

incoming boundary conditions must be selected for the corresponding source, Y m (u). We 
therefore substitute Y n (u) = (1 — -u 2 )~ jro / 4 F m (u) where F m (u) is regular at the horizon. Since 
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F m (u) is regular at the horizon we may expand in a power series in to and solve order by 
order 4 



Y n {u) = (1-u 2 )-^ x (4.18) 
1 + — {tt - 4tan- 1 ( v ^I) - 6 log 2 + 41og(l + y/u) + 21og(l + u)} + 0(m 5 



The insertion of this solution into the action reduces to a boundary term. After re-inserting 
units ^ 

y u (u) = (nT) 2 y n (tv = —) Y w (uj, u) = Y n (rv = — , u) , 

and denoting A(u) as the coefficient in front of the kinetic term as in Refs. jl^,!^] we have 

R 2 (ttT) 2 f 

G R {u) = -A(^y_ w (i03„y w (iOUo, A(u) = -^-p- — ye ■ ( 4 - 19 ) 

Substituting the solution Eq. ()4.18|) into this expression we obtain 

— 2T 

k = lim lmG R (u) , (4.20) 

= V\T 3 n. (4.21) 

The divergence which arises as u — > does not effect the imaginary part. In the last step 
we have used the relation R 2 /a' = \^X. The significance of the diffusion computation will 
be discussed in the next section. 



V. DISCUSSION 

In summary we have computed the heavy quark diffusion coefficient in M = 4 SYM 
by exploiting the AdS/CFT correspondence. From Eq. (J4.20)) and Eq. (|2.Hj) the diffusion 
coefficient is 

D = ^ (5.1) 

= ^7T (5 ' 2) 

where A = g\ M N c . This result is interesting both theoretically and phenomenologically. 

An immediate observation is that the heavy quark diffusion coefficient D is parametri- 
cally small compared to the momentum diffusion coefficient r)/(e+p) = l/(47rT) [23] , i-e. D 
depends on the peculiar, but characteristic, l/y/X. This should be contrasted with pertur- 
bation theory where all transport scales are the same order of magnitude, 1/(A 2 T) log(A _1 ). 

Theoretically, the diffusion computation proceeded as follows. A Wilson line running 
along the Schwinger-Keldysh contour is the partition function of a single heavy quark 5 . 



4 We thank A. Starinets for providing us this expression. 

5 The shape of the contour is irrelevant. Taking a contour which runs straight down the imaginary time 
axis to — i/3 gives the Polyakov loop. This is the more common but equivalent definition of the heavy 
quark partition function. 
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There is a Wilson line on the "1" axis and a Wilson line on the "2" axis. In a Minkowski 
metric the corresponding semi-classical string solution must span the full Kruskal plane if 
the right (left) quadrant of the Kruskal diagram is to be identified with "1" ("2") fields in 
the real time path integral. We find this semi-classical solution in Section II VI and it is a 
single straight string in AdSc, x S$ which extends from the boundary to the event horizon 
when viewed from an observer in the right quadrant. Force- force correlators in the presence 
of a heavy quark are electric field correlators with links running along the contour after the 
heavy quark fields have been integrated out. These same correlators may be obtained by 
variations of the unperturbed Wilson loop. Thus by varying the end point of the relevant 
string solution in AdS x S§ one may compute real time force-force correlators, i.e. the 
diffusion coefficient. 

Given the diffusion coefficient in Eq. (j5.1|) . the quark must be sufficiently heavy for the 
Langevin theory to apply. Since the Langevin theory is classical, consistency demands that 
the relaxation time be large compared to the inverse temperature 



%D»L. (5.3) 



This leads to the following constraint 



M^>— VX. (5.4) 

For a string of length L stretching outward from the horizon we have, M ~ L/a', 
ttT = r /R AdS , and y/X = R 2 AdS /a'. This Langevin constraint is satisfied by the gravity 
computation whenever 

L>r . (5.5) 

i.e. whenever the length of the string is long enough. Clearly mass or energy of order Ty/X 
sets a boundary for a quasi-particle picture to be valid. In the gravity computation this 
corresponds to a string of length L ~ r . 

Phenomenologically the diffusion coefficient is an extremely interesting number. Bearing 
in mind that TV = 4 SYM is not QCD, the authors still want to substitute numbers 

0-9 / 1.5 y /2 

2ttT \asYMNj ' 1 • ' 



This could be compared to extrapolations of weak coupling to QCD jy,|4 



6 / 1.5 x 2 



2ttT \a s N c 



Of equal importance is the mass scale where we expect the heavy quark theory to apply 

1/2 



(a^Gev) 



(5.8) 



This suggests that the Langevin process might not be applicable to charm quarks, though 
of course there is an unknown proportionality factor in this equation. 
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The RHIC data favor a diffusion coefficient of about 3/(2ttT). However, relativistic 
effects obscure the relation between the diffusion coefficient and the semi-leptonic data. The 
role of radiative and collisional energy loss remains unclear 0, ^j. The next step from the 
perspective of the AdS/ CFT correspondence is to consider quarks with finite velocity - work 
is in progress in this direction [44]. The framework set up here should allow a computation 
of the real time transport properties of these quarks. 

Note Added. In the few days surrounding this preprint three similar studies of heavy 
quark energy loss appeared. The first of these computed the transverse momentum diffusion 
of an ultra- relativistic quark j4^|. The second paper, denoted HKKKY after the authors, 
computed the drag on a quark moving with finite velocity |4f|. A third paper j4?| indepen- 
dently computed the drag using methods similar to HKKKY and obtained the same answer. 
The setup in these papers differs signficantly from this work. 

The HKKKY paper is particularly instructive. Consider the Langevin process of a heavy 
quark 

^ = m - VDPi , mM*)) = «w - o . ( 5 - 9 ) 

where the drag, fluctuations, and diffusion are related by Einstein relations 

k 2T 2 
VD= 2MT> D = ~- 

In this paper we considered times short compared to the relaxation time (M/T) D, and 
computed the strength of the noise k, taking the mass to infinity. The mass is not needed 
to compute the diffusion coefficient or dE/dx = T/D. 

HKKKY computed the drag by considering the change in the average momentum over a 
time longcompared to the time scale of noise correlations, ~ 1/T. The Langevin process 
(see e.g. p) predicts the probability that a quark with momentum p at time will arrive 
with momentum p at later time t 

-P(p ? *IpojO) = — exp 
V ' ; Zj2nMT (1 - e~ 2 ^) 

The average momentum obeys 

(p(t)> = poe-to* , (5.10) 

while the width obeys 

<(Ap) 2 ) = 3MT (1 - e~^ Dt ) . (5.11) 

HKKKY considered an atypical quark, with momentum p ^> V MT, and computed the 
change in the average momentum over a time interval of order, 1/rjjj = (M/T)D. This 
determines the drag, rj£>. For an atypical quark the fluctuations are neglible over this time 
interval. They also compute the mass and then are able to deduce diffusion coefficient 
through the Einstein relations. The result agrees with Eq. (|5.1|) : the AdS /CFT corre- 
spondence is neatly consistent with the fluctuation dissipation theorem. It would be quite 
interesting to see the full structure of the Langevin Green function emerge from the string 
theory. In essence this computation has been performed already through an amalgamation 
of our works. 



(p - ppe-^y 

2MT(l-e" 2 ^) 
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HKKKY also computed the heavy quark energy loss at finite velocity. The bending string 
solution they obtained agrees with our preliminary computations. 

Acknowledgments. We wish to thank Andrei Starinets and Edward Shuryak for fruitful 
discussions. JC thanks M. Kulaxizi, S. Giombi, R. Ricci and D. Trancanelli. We especially 
thank Ismail Zahed for challenging us to find a clean formulation of the problem. This work 
was supported by grants from the U.S. Department of Energy, DE-FG02-88ER40388 and 
DE-FG03-97ER4014. 



APPENDIX A: THE HEAVY QUARK PARTITION FUNCTION 

In this section we clarify several aspects concerning the heavy quark partition function. 
We follow closely the discussion of McLerran and Svetitsky [36| • We introduce the operators 
Q(x, t), Q^(x,t) which create and annihilate static quarks at point x and time t. These 
fields satisfy the anti-commutation relations 

{Q i (x,t),Qt( x ',t)} = <^(x-x') . (Al) 

From the Lagrangian Eq. (j3.1J) the time evolution of quark fields is given by 

Q(x, t) = e~ lM ^U(t, * )Q(x, t ) , (A2) 

where both t and to axe in the Schwinger-Keldysh contour and U(t,to) is the transporter 
between these two points at fixed position x. As in Ref. j^, the partition of function in 
the presence of a heavy quark is 



J2( s \ e ' PH \ s ) = [ d 3 xJ2(s'\Q^-T)e-P H QHx,-T)\s>) , (A3) 

s s' 

= / d 3 xJ2(s' |e-^Q(x, -T - iP)Q\x, -T)\ s') , (A4) 

J a' 

in which \s) is a state of the system with only one heavy quark, and \s') is a state with no 
heavy quarks (i.e. Q\s') = 0). By means of Eq. (jA2|) we obtain Eq. (|3.18j) . With this 
definition of the partition function, the force-force correlator is 

D>(t) = (mm) HQ = ^-J2( s \ e ~ PH:F (t)Hu)\s) ■ (as) 

We introduce a complete sets of states between the two force operators. We notice that, 
since we are assuming weak fields (which cannot create heavy quarks) and since the force 
operator does not change the number of quarks, it is enough to consider a complete set of 
one particle heavy quark states plus bath. Following Ref. [34| we obtain the KMS relation, 

£>>(*) = (F(0)F(t + i/3)) HQ = D<(t + t(3) . (A6) 

Since the force is a Hermitian operator local in time, the definition of the spectral density 
pit) = D > {t) — D < (t) leads to the standard relations between correlators, Eq. (j3.11j) . (j3.12|) 
and (I3~T3J) . 
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We conclude by remarking that, since the states \s') can be considered as the |0a) Fock 
states for the heavy quarks in a gauge field background, the partition function is 

Z HQ = J d'x Tr [(0 A \ Q(x, —T - i/3)Q\x, —T) |0 A >] , (A7) 

where the Tr represents the thermal trace over the bath (the Yang-Mills fields). Thus, 
representing the vacuum expectation value in Eq. (jA7|) as a path integral we arrive at 
expression Eq. (J3.15)) . The effect of the vacuum average is to introduce the ie prescription 
in Eq. 1)3.16)1 as in zero temperature field theory [isl ]. 
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